o 

(N 
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The quantum Hall transition [l|] is one of the simplest and most studied quantum phase transitions. 
Nevertheless, the experimental observation of a new phase in this regime, the quantum Hall insulator, 
still remains a puzzle since the first report more than a decade ago 043 1 ^^ i* is in contradiction 
with all theoretical studies based on microscopically coherent quantum calculations [9l-[l2]. In this 
work we introduce into the coherent quantum theory a new ingredient - rare incoherent events, 
in a controlled manner. Using both direct numerical solutions and real-space renormalization, 
we demonstrate that these decoherence events stabilize the elusive quantum Hall insulator phase, 
which becomes even more stable with increasing temperature and voltage bias, in agreement with 
experiments. 

PACS numbers: 73.43.Cd,73.40.Hm,72.20.My 
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The integer quantum Hall effect [l| has been a 
paradigm for two-dimensional quantum phase transi- 
tions: a transition between the quantum Hall phase, 
characterized by a quantized Hall resistance pxy and a 
vanishing longitudinal resistance pxx, and an insulator, 
characterized by diverging pxx and pxy- This transition 
can be intuitively understood within the semiclassical 
description of the quantum Hall effect, valid at strong 
magnetic fields. In such fields electrons follow equipo- 
tential lines, which due to disorder arc localized around 
valleys of the potential, at low energies, and peaks, at 
high energies. Near the potential saddle points, where 
such trajectories get close to each other (Fig. [T^), elec- 
trons can tunnel from one trajectory to another, where 
the tunneling probability depends on the characteristics 
of the saddle point [l3|. The critical energy is the en- 
ergy at which there will be a trajectory that percolates 
through the system. Thus the quantum Hall transition 
may be described as a quantum percolation transition 
[14l4l6| . Such a network model [17[ of saddle points has 
also formed the basis for extensive numerical calculations 
and for real-space renormalization group (RSRG) calcu- 
lations [18[. These calculations, consistent with other nu- 
merical studies of the quantum Hall transition [19[ and 
with field theoretical renormalization group analysis [20[ 
indeed demonstrate the existence of a critical value Tc 
of the average transmission through a saddle point T, 
separating two stable phases - the quantum Hall phase 
[T > Tc) and the insulating phase (T < Tc)- 

However, the experimental discovery of a new phase, 
the quantum Hall insulator [2-|8[, characterized by a 
quantized pxy but an exponentially large pxx (inset of 
Fig. 11^). has eluded theoretical understanding. Micro- 
scopically coherent quantum mechanical calculations, ap- 
plied to this regime, such as numerical simulations [10[ 
or RSRG analysis [IjJ, have demonstrated that there is 
no quantum Hall insulator phase in this system, consis- 
tent with other studies [9|, ll^l , except perhaps for small 
enough systems J2l| . Taking dephasing into account phe- 
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FIG. 1: (color online) (a) Electron trajectories in the po- 
tential landscape. In strong magnetic field electrons follow 
equipotential lines (thick lines), and may undergo incoher- 
ent scattering events (wiggly lines). Near saddle points an 
electron can tunnel from one equipotential line to another 
(broken curves), with a corresponding tunneling probability 
T. Such a junction, including the possibility of quantum tun- 
neling and of incoherent scattering, is represented by the el- 
ementary unit (b), characterized by T and the decoherence 
parameter e. These elementary units, represented by orange 
dots, are connected in a network (here of size 3 x 3). (c) The 
renormalization group procedure: mapping five elementary 
units connected in Whetstone bridge form, representing the 
system in (a), onto a single effective elementary unit. 



nomenologically, by addressing a finite-size system, mim- 
icking the finite dephasing length, Pryadko and Auerbach 
[lOj have concluded that pxy of the entire system will be 
given by its value for that length scale, which, in princi- 
ple, could be exponentially large - again in contradiction 
with the experimental observation. 

Here we incorporate into the quantum mechanical cal- 
culation rare incoherent scattering events the electron un- 
dergoes during its motion along the equipotential lines 
(wiggly lines in Fig. [1^). These events are introduced 
into the model via phase randomizing, current conserv- 
(see Fig. Wp), which means that for 
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mg reservoirs 

every electron that enters such a reservoir there is an 
electron that leaves it. However, the phase of the out- 
going electron has no correlation with the phase of the 
incoming electron, so that once an electron enters such 
a reservoir, interference effects are destroyed. The basic 



unit in our network model consists of a saddle point (de- 
scribed by a scattering matrix, with random phase due to 
disorder) straddled by two current-conserving reservoirs, 
with probability e to enter each reservoir (Fig. [T|d). Ac- 
cordingly, the transport parameters of the system can be 
calculated via a standard quantum multi-terminal scat- 
tering approach |23| . where two of the terminals repre- 
sent the incoming and outgoing current, two - the volt- 
age probes which measure the Hall voltage, and the rest - 
the current-conserving reservoirs, leading to decoherence. 
For e = this model reduces to the previously studied 
coherent network model 171. The coherent transmission 



from left to right, Tr^, is thus the quantum probability 
that the electron traverse the system in this direction, 
without entering any of the phase randomizing reservoirs. 
Tlr, Rl (the reflection from left to left) and Ru are simi- 
larly defined. Consequently, the probability e to enter the 
reservoir determines the amount of decoherence in the 
system, characterized, e.g. by the decoherence length L^, 
the typical length an electron travels before losing mem- 
ory of its original phase. When e = 0, no electron enters 
the phase randomizing reservoirs, and the system is fully 
coherent, L^ — )> oo. When it is unity the system is classi- 
cal, with no coherent transport. Tuning this probability 
from zero to unity allows us to probe the crossover be- 
tween the fully coherent to the fully incoherent regimes. 

We first solve numerically for p^x and pxy for a net- 
work of size L X L (Fig. [T}3) , for different values of the 
average transmission though a saddle point, T and the 
decoherence parameter e. pxx = (1 — T)/T is determined 
by the effective transmission T through the whole sys- 
tem, while Pxy, as in the experiment, is determined by 
Vh, the anti-symmetric component of the difference be- 
tween the chemical potentials at the upper and lower 
branches of the structure with respect to magnetic field, 
which is nonzero due to the chiral nature of the problem: 
Pxy = Vh/I, where / is the current. Because both pxx 
and Pxy are exponentially distributed, we have used a log- 
arithmic average [ll| to calculate the effective renormal- 
ized values. The values of the saddle-point transmission 
probabilities T are taken from a wide distribution, with 
a predefined average, while the coupling to the current 
conserving reservoirs e are taken from a delta distribu- 
tion. 

In Fig.[2]wc plot pxx (a) and pxy (b) as a function of the 
size of the system, for different values of the dephasing 
parameter e. In the insulating phase (T < Tc), and in the 
absence of decoherence events (e -^- 0), both pxx and pxy 
increase exponentially, consistent with previous studies 
[ld |. In the presence of decoherence, pxx first increases 
with system size (for L < L^), and then saturates, as 
one expects for a classical system. Surprisingly, while 
Pxy also initially increases, for L < L^, it reaches a max- 
imum and then decreases. For the samples with larger e 
(smaller L^) pxy decreases all the way to unity (all re- 
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1. (d) Results of the real-space renor- 
Pxy (inset Pxx) as a function of system 



FIG. 2: (color online) pxx (a) and pxy (b) from the network 
model as a function of network size for several values of the 
dephasing parameter e (values depicted in (b)), and average 
transmission of a single scatterer < T >= 0.13. While pxx 
saturates at L 3> 1/^4, Pxy exhibits a nonmonotonic behavior. 
(c) Same date as in (b), plotted as pxy — 1, multiplied by a 
constant, as a function of size on a double-log scale. At large 
sizes all curves collapse onto a single curve, confirming that 
for large sizes pxy 
malization method 
size for different values of the initial dephasing parameter e , 
with the same initial transmission < T >= 0.4, exhibiting the 
same behavior of both pxy and pxx as the network model. In 
the limit L ^ L^, pxy is quantized to unity, while pxx can be 
arbitrarily large. This is the quantum Hall insulator phase. 



sistance values are expressed in units of /i/e^, where h is 
the Planck constant and e the electron charge). While 
for samples with smaller e, pxy has not yet reached the 
asymptotic regime, L ^ L^, we demonstrate in Fig. [2j; 
the collapse of all the curves, when plotted as pxy — 1 vs 
L. This indeed confirms that, independent of e (or L^), 
Pxy scales as pxy = l + c^f (L/L^), with f{x) ~ l/i/i for 
large x. This phase, where pxy is quantized to unity, and 
Pxx could be exponentially large, is the elusive quantum 
Hall insulator phase. 

In order to address asymptotically large systems L 3> 
i^, we employ the RSRG method. In this approach one 
replaces a part of the system, containing several saddle 
points, by one effective saddle point, whose character- 
istics depend on those of the saddle points included in 
that subsystem (Fig. ^). By following the dependence 
of the effective saddle-point transmission probability, T, 
upon renormalization (see below), one can demonstrate 
that the coherent system flows either toward the insu- 
lating state (T = 0) or towards the quantum Hall state 
(T = 1) [ll|, an observation consistent with field theo- 
retical renormalization group analysis [20| and numerical 
calculations [19[. In the presence of decoherence, both T 
and e are renormalized, as the probability for an incoher- 
ent scattering event increases as the system size increases. 
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FIG. 3: (color online) (a) p^y and pxx for ordered system 
(same phase and transmission value for all the scatterers) as 
a function of length for system of fixed width W — 4, the 
green color is for T = 0.05 and the blue is for T — 0.1. Here 
dephasing occurs only at the corners of the elementary squares 
(4x4 networks, plotted as empty dots) for which e = 1. p^x 
grows linearly with the length while pxy saturates, (b) For a 
system of fixed length, as the system becomes wider pxy goes 
to unity {Cxy equals 43 and 17.6 for T = 0.05 and T = 0.1 
respectively), pxx goes to zero as expected, with Cxx equals 
3500 and 718 for T = 0.05 and T = 0.1, respectively. 



The basic cell for the RSRG method used here is a five- 
unit cell in a Wheatstone bridge setup (representing the 
system depicted in Fig. [1^) , which is transformed into a 
single unit virith effective transmission Teff and decoher- 
ence parameter Ce//- Ce// is defined in terms of all the 
coherent parts of the scattering matrix. 



Ee// 



= 1 - 



Tlr + TjiL + Rr + Rl 



(1) 



Our RSRG transformation thus consists of two recursion 
equations, one for the total current (transmission) and 
one for the decoherence parameter e. We solved these 
equations using Monte-Carlo sampling rejection method 



31| . and averaged over 2.5 * 10^ realizations. Starting 
from some initial distributions of T and e, this step is 
repeated to generate the next generation distributions, 
which are used as input to the next iteration. Fig. [2]i 
depicts the resulting pxx and pxy as a function of the 
system size L, for a given value of T, and for several 
values of L^ (or e) on the insulating side. Similarly to 
the direct numerical solution, p^x saturates for L > L^, 
and Pxy displays the same intriguing behavior: initially it 
grows exponentially, but then around L ~ L^ it exhibits 
a maximum and then decreases back to unity. 

While the decrease of pxy towards unity, as the system 
size increases, seems apriori surprising, one can show that 
it is a direct consequence of the rules of connecting re- 
sistors in series and parallel. To demonstrate this point 
we have calculated Rxx and Rxy (the longitudinal and 
Hall resistances for non-square systems) for a stack of 
coherent ordered squares, each of size LqX Lo, where the 
only decoherence scatterers are at the corners of these 
elementary squares (see inset of Fig.[3|), connected either 
in series (Fig. [3K) or in parallel (Fig. [Sb). For the series 
connection we find that Rxx increases linearly with the 
system length L, as expected, and Rxy saturates, while 



for the parallel connection, both Rxx and Rxy decreases 
as 1/W, where W is the width of the sample, the former 
towards zero, again as expected, while the latter towards 
unity. This same behavior is also observed for rectangu- 
lar disordered network- model systems (not shown). Both 
these behaviors of Rxy can be readily understood. We 
first note that since / = TV, where V is the voltage dif- 
ference between source and drain, and Rxx = (1 — T)/T, 
one can write i?xj/ = {VH/V){Rxx + i)- For the series con- 
nection, when L ^ L^, one can think of the system as 
consisting oiL/L^ coherent segments, connected incoher- 
ently. Thus the voltage drop on each segment is VL^/L. 
Because the Hall voltage of each segment is linearly de- 
pendent on the voltage drop across that segment, it scales 
like 1/ L, and as Rxx grows linearly with L, Rxy remains 
constant. On the other hand as the system width W in- 
creases, for constant length, Rxx decreases as 1/W, and 
thus the above relation dictates that Rxy also decreases 
(since Vh is bound from above by V). In fact, in this 
limit, the upper chemical potential becomes dominated 
by the source chemical potential, while the lower chemical 
is dominated by the drain chemical potential, and thus 
Vh/V approaches unity as the width increases. Conse- 
quently, as Rxx approaches zero with increasing width, 
Rxy approaches unity. (This is in contrast with the anal- 
ysis of Ref. lid , which claims that Rxy remains constant 
as the width increase, while Rxx decreases towards zero, 
violating the above relation between Rxx and Rxy) Pxy 
for large square of size L x L can then be obtained by 
first making the system longer, of length L, such that its 
Rxy does not change any more, and then increasing its 
width to L, so that the Hall resistance decreases towards 
unity. This is a generalization of the case considered in 



Ref. I24l which effectively connected in series and paral- 



lel puddles of Rxy = 1, leading to pxy — 1 for the full 
system. This observation is also consistent with the two- 
phase approach [25| . 

Having established the stability of the quantum Hall 
insulator phase, we compare our results, using both 
methods, in Fig. 31 to the experimental data. Panels 
(a) and (b) depict our calculation and the experimen- 
tal data (inset of (a)), pxx and pxy are plotted as a 
function of the distance from the critical point (T — Tc 
in the theoretical curves, B — Be in the experimental 
curve). Several theoretical curves, for a system of fixed 
L, but of different e (or L^) arc plotted, demonstrating 
that the quantization of pxy in the quantum Hall insu- 
lator phase becomes more exact as the level of decoher- 
ence increases (larger e, smaller L^). Interestingly, the 
experimental curves exhibit better quantization with in- 
creasing temperature, which we attribute to increased 
decoherence. Enhanced decoherence also explains the 
better quantization of pxy for higher currents [j|. For 
even higher temperatures, approaching the energy gap 
in the quantum Hall regime, one observes breakdown of 
the quantization in both phases - the quantum Hall phase 




FIG. 4: (color online) Theoretical (a,b) and experimental [21 
(inset of (a)) results for the longitudinal and Hall resistivity 
plotted as function of the deviation from the critical point 
for systems with different phase lengths (theory) or tempera- 
tures (experiments). The quantization of the Hall resistance 
on the insulating side improves with increasing decoherence 
(theory) or temperature (experiment). Symmetry between 
the quantum Hall phase and the insulating phase, p^x on the 
insulating side and 1/ Pxx on the quantum Hall side, experi- 
ment [31 (inset of (c)) and theory (c,d). Theoretical results 
are for system with finite phase length. The data is plotted 
with the axis for 1/pxx inverted (top axis), for two differ- 
ent temperatures (experiment) and phase lengths (theory), 
demonstrating the symmetry between the two phases. The 
inset demonstrates that this symmetry is obeyed in the the- 
ory for many orders of magnitude. 



and in the quantum Hall insulator phase [J-ly, |8|. An- 
other striking feature of the experimental data [3 was the 
symmetry of pxx on the two sides of the critical point, 
Pxxi^v) = 1/ pxx{—^v), where v is the filling factor, the 
number of electrons in the system per available states in a 
Landau level, inset of Fig.|4l:. This symmetry is also man- 
ifested in our results (Fig. [3J:,d). In the coherent case it 
can be traced to the symmetry of the disorder potential, 
which leads to7^ = l — 7^ = 1/2. In that case it is easy 
to see [12] that since, by definition, in the fully coherent 
case pxx{T) = T/(l - T), then clearly pxx{Tc + AT) = 
(Tc + Ar)/(1 -Tc- AT) = l/pxx{Tc - AT). In the 
presence of incoherent scattering, pxx is given, as dis- 
cussed above, by the coherent p^x on a scale of L^. Since 
in the model the decoherence parameter is, by definition 
(Eq. [1]) symmetric around the critical point, then the 
above relation is obeyed perfectly (see inset of Fig. [4ji). 
The experimental deviations from this symmetry (inset 
of Fig. [4|;) thus makes it possible to investigate the de- 
pendence of the decoherence length on magnetic field and 
density, allowing a deeper understanding of the nature of 
the incoherent processes at such low temperatures. 
The relevance of incoherent scattering at milli-kelvin 



temperatures in the quantum Hall regime, imperative, 
as shown here, for explaining the quantum Hall insulator 
phase, has already been established experimentally [2g . 
Incoherent scattering should be explored in the context 
of other quantum phase transition as well. In particu- 
lar, it may also explain other puzzling two-dimensional 
phenomena, such as the apparent metal-insulator transi- 
or the intermediate metallic phases observed 
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tion 

in the superconductor-insulator transition in disordered 
thin films [29j , and in the quantum Hall to insulator tran- 
sition [30|. The present calculation allows quantitative 
determination of the incoherence length which is crucial, 
for example, for any possible application of mesoscopic 
devices as quantum bits, the basic building blocks of a 
quantum computer. 

We thank A. Auerbach and A. Stern for fruitful dis- 
cussions. This work was supported by the ISF and BSF. 
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